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ABSTRACT 

We perform a self-consistent calculation of the thermal structure in the crust of a superbursting neutron star. 
In particular, we follow the nucleosynthetic evolution of an accreted fluid element from its deposition into the 
atmosphere down to a depth where the electron Fermi energy is 20 MeV. We include temperature-dependent 
continuum electron capture rates and realistic sources of heat loss by thermal neutrino emission from the 
crust and core. We show that, in contrast to previous calculations, electron captures to excited states and 
subsequent 7-emission significantly reduce the local heat loss due to weak-interaction neutrinos. Depending on 
the initial composition these reactions release up to a factor of 10 times more heat at densities < lO" gcm^^ 
than obtained previously. This heating reduces the ignition depth of superbursts. In particular, it reduces the 
discrepancy noted by Gumming et al. between the temperatures needed for unstable '^C ignition on timescales 
consistent with observations and the reduction in crust temperature from Cooper pair neutrino emission. 
Subject headings: dense matter — nuclear reactions, nucleosynthesis, abundances — stars: neutron — ^X-rays: 
binaries — X-rays: bursts 



1. INTRODUCTION 

The ability to regularly monitor the X-ray sky with instru- 
ments such as RXTE and BeppoSAX has produced many ex- 
citing discoveries in the last few years. An excellent example 
are superbursts (for a review see KuuUcers 20 041). Like normal 
type I X-ray bursts (for a review, see iStrohmaver & BildstenI 
|2004), superbursts are characterized by a rapid rise in the 
light curve followed by a quasi-exponential decay; when com- 
pared to type I X-ray bursts, superbursts are roughly ^10^ 
times more energetic, have cooling timescales of hours, and 
recur on timescales of years. For one source, 4U 1636—53, 
three s uperbursts have been observed over a span of 4.7 yr 
dWiina nds 2001; lKuulkers et al.ll2004l) . Although most super- 
bursts are observed at mass accretion rates 0. 1-0.3 Mfidd, with 
MpAd ~ 10^^ gs~' bei ng the Eddington accretion rate, recently 
lin 'tZand et al.l ( l2004 a) detected superbursts from the rapidly 
accreting source GX 17+2. 

The currently favored scenario for these superbursts 
is the thermally unstable ignition of ^^C at densi- 
ties ~ 10^-10^ ^. ' r"^^. (Gumming & BildstenI 
IStrohmave r & Brownl l2002h . 
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ICumming & MacbethI (120041) 
demonstrated that the superburst lightcurve after the peak 
luminosity has been reached declines as a broken power law, 
with the early behavior dependent on the burst energy and the 
later behavior depen dent on the igniti on depth. This behavior 
has been observed (ICumming et al.l 120061) . The ability to 
constrain the burst energetics and ignition depth is crucial, not 
only for understanding the ignition of ^^C in a heavy-element 
bath, but also for p robing the inte r ior of t he neutron star 
ferown 2004; Co oper & NaravanI 120051: ICumming eta[\ 
I2OO61) . iBrowni (2004) demonstrated that the temperature at 
densities ~ 10^ gcm^^ is sensitive to the thermal conductivity 
of the inner crust and the neutrino emissivity of the core; 
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this was confirmed by ICooper & Naravan' ('2005') in a more 
extensive study. Recently, [Cumming et al. (2006) showed 
that the neutrino emission from singlet-state Cooper pairing 
of neutrons in the inner crust limited the temperature at 
neutron drip to be < 5 x 10^ K. Taken at face value, this 
would lead to superburst ignition column depths about an 
order of magnitude deeper than the 0.5 — 3 x lO'^ gcm^^ 
deduced from observations, and correspondingly to recur- 
rence times that ar e about an order of magn itude longer than 
observed. This led lPage & Gumming! (l2005h to speculate that 
the superbursters might in fact be strange stars. 

The unstable ignition of '^C depends critically on the ther- 
mal structure of the outer crust and is sensitive to the heat 
sources located there. The outer crust of the neutron star 
is composed of nuclei and degenerate electrons and exists 
where the electron chemi cal potential ji^ < 30 MeV (see 
iPethick & RavenhaiHl 19981 for a lucid discussion). The ashes 
of hydrogen, helium, and carbon burning in the neutron star 
envelope are continuously compressed by the ongoing accre- 
tion of matter and incorporated into the crust. The increas- 
ing density induces nuclear reactions at various depths that 
can release energy. We show that most of this energy is de- 
posited in the crust, rather than being carried away by neu- 
trinos. This leads to a hotter crust and decreases superburst 
ignition depths and recurrence times. 

Previous studies of the crust heat sources used a single 
representative i sotope, ^^Fe dSatol 119791; iHaensel & ZdunikI 
119901) or '"'^Pd (Haensel & Zdunik '2003') for the product of 
H/He burning. Haensel & Zdunik (1990) found that elec- 
tron captures would occur in two stages (see Fig. [TJ. The 
first electron capture (labeled "1" in Fig. [T|i would occur 

"tiw o.s-gs capture onto an even-even nucleus: 
(Z — 1 , A) + v^., where Z and A are even. Here 
(c, where the threshold for 



when pic 

(Z,A)+, 

the threshold is /Sthr = ^thr.gs-gs 
the ground-state-to-ground-state transition /ithr.gs-gs is com- 
puted from atomic mass differences and therefore includes the 
electron rest mass, and where Eexc is the energy of the lowest 
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excited state into which the capture can proceed. Because of 
the odd-even staggering of the nuclear masses, E^^^ ^^^^ for 
electron capture onto the resulting odd-odd nucleus (Z — 1 , A) 



is less than E, 



(ZA) 



> '^tto^glgi ^"'l ^ sec- 



-thr,gs-gs- Consequently, 
ond electron capture (Z — 1,A) + ^ (Z — 2, A) + Ve im- 
mediately follows. Where experimental data were lacking, 
which is the case for the vast majority of nuclei, the heat de- 
posited into the crust Q was calculated assuming that the cap- 
ture was into the ground state of the daughter nucleus (Fig.[Tl 
lab el "4"). The heat deposite d into the crust is then estimated 
as dHaensel & ZdunikI 120031) Q = (/ig - £'thr,gs-gs)/4, where 
the factor 1/4 arises from the integral over the reaction phase 
space in the limit that (/Xf — -Ethi) //X^. <C 1 . This estimate as- 
sumes that final states are very low in excitation energy so that 
the contribution to the heat deposition from the radiative de- 
excitation of the daughter nucleus is always negligible com- 
pared to the heat release from the electron capture. In this 
picture, 75% of the nuclear energy release is emitted in form 
of neutrinos. 
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Fig. 1 . — Electron capture for which the first transition (labeled "1") pro- 
ceeds to the ground state of nucleus ''Z — 1, but the second transition instead 
goes to an excited state "2" followed by a radiative de-excitation "3". Previ- 
ous works assumed that the second transition always went to the ground state 
"4". 



In this paper, we calculate heating from nuclear reactions 
in the outer crust (where the mass density is less than neutron 
drip, ^ lO'^ gcm^^) using realistic electron capture rates ob- 
tained by microscopic calculations of the transition matrix el- 
ements between the parent ground state and daughter excited 
states multiplied with the phase space factor corresponding 
to the Q-value of the decay. The rates depend on the exci- 
tation energy and the ambient temperature and density. We 
use a full reaction network having temperature- and density- 
dependences. Our thermal model of the crust, described in 
§ 12] is coupled to the reaction network (§ |3]l, so that our so- 
lution is self-consistent. Instead of starting with a single nu- 
cleus, we follow the actual mix of nuclei initially produced 
by surface burning, both from the rp-process (§ 13.2.1b and 
a mixture in nuclear statistical equilibrium formed during a 
superburst (§ |3.2.2| l. We find that electron captures deposit 
considerably more energy in the crust than in previous esti- 
mates (Ha ensel & Zdunik ^003). This larger heat deposition 
is mainly due to the electron captures populating excited states 
with energy iSexc (Fig.[T] label "2"), rather than capturing into 
the ground state as assumed previously. The populated ex- 
cited states then radiatively de-excite (Fig. [T] label "3"), de- 
positing the entire excitation energy as heat into the crust. We 
also introduce an approximate model that gives a good esti- 
mate of the heating in the outer regions of the crust for an 



arbitrary mixture in § 13.31 and describe the impUcations for 
superburst ignition in § |4] 

2. THE THERMAL STRUCTURE OF THE CRUST 

Our thermal model is the same as used by iBrownl (l2000i 
l2004l) and we shall only summarize its essential proper- 
ties here. Over the crust, the gravitational potential <I> « 
jc^ln(l — 2GM/rc^) is nearly constant; moreover, we are al- 
ways in hydrostatic equilibrium and can approximate the pres- 
sure as being nearly independent of temperature. Under these 
approximations, we then solve the equations for temperature, 
luminosity, and nuclear heating. 



L(l+z) 



dr 4Kr^K 

_2^/,2 d(g2'^A-'L) dLnuc 

' dr 
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AnP-{\ +z)ev, 
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(2) 
(3) 



where z = (1 — 2GM/rc^) '/^ — 1 is the gravitational redshift, 
Q{r) is the heat deposited per accreted nucleon as computed 
from our reaction network (§[3]l, and nia = 1.66 x 10^^^ g = 
1 u. Over the crust the potential <I> varies only slightly and our 

code simplifies equations ([T]i and (|2]i by holding e*/' con- 
stant. 

We model the core with an analytical density pre- 
scription p = po[l — (r//?)^] dTolm an 1939), where po — 



\-2 



-1 



. Here 



15(87r)-'M/?-3 with R = {2GM/c^) 1 - (1 

r and R are respectively the circumferential radial coordinate 
and stellar radius. This prescription gives a reasonable ap- 
proximation to the p (r) that would be obtai ned from a mod- 
ern n uclear equation of state (EOS) (see iLattimer & PrakasE 
1200 iL and references therein). With this choice for the 
EOS, the gravitational mass of the neutron star M and red- 
shift z specify the mechanical structure of the star. We 
choose the crust-core boundary to be at a transition den- 
sity 1.6 X lO'^'^gcm^^, which is the density found by a 
Max well construction us ing an EOS for the inner crust 
(Nege le & VautherinlfT9T3 ) and the core EOS jAkmal et alj 
il998i) . We then adopt a neutron star of mass 1.6 M©, and 
set the radius at the crust-core boundary to 10.5 km. This 
choice makes [1 — 2GM/(rc^)]^'''^ = 1.35 at the boundary. 
Our neutron star then has a radius 10.8 km and a surface grav- 
ity 2.43 X lO'^^cms-^. 

For our numerical calculations we use an accretion rate, 
measured in the rest frame at the surface, of M = 3.0 x 



; the local accretion rate, per unit area, is then 2.1 x 

cm^^ 



10"gs-' 

10^ gs^' cm^^ and the luminosity (for spherical accretion) is 
0.3 of the Eddington luminosity for a solar composition^. The 
core neutrino emissivity is modified Urea with an emissivity 
10^''(r/10'^ K)** ergscm^^ s^'. Note that because our starting 
density for the integration is p = 6.2 x 10^ gcm^^, which is 
just below where the rp-process burning has concluded, we 
do not resolve the H/He burning shell in the envelope. Ob- 
servations of the ratio of persistent fluence to burst fluence 
sugge st that a good fract ion of the accreted H/He burns sta- 
bly din 't Zand et alj |2004b). so we set the temperature at this 
boundary to be 4.2 x 10** K, which is computed by integrating 

^ Here M is the baryon accretion rate expressed in mass units, and not the 
rate at which the gravitational mass of the star changes 
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steady H burning dSchatz et al.ll 19991) at a local accretion rate 
ofm = 2.6 X 10 gs^' cm^^. 

Solving equations ([T]l and (|2]l requires the EOS, thermal 
conductivity K a nd neutrino emiss ivity gy. We use a tabulated 
electron EOS (Timmes & Swesty 2000) with ion electrostatic 
interactions (Farould & Hamaguchi 1993i). The pressure of 
the free neutrons is computed using a zero-teniperatu re com- 
pressible liquid-drop model jMackie & Bavmll 19771) . Heat 
is transported by degenerate, relativistic electro ns; our treat- 
ment of the electron scattering frequency follows' Brown et al] 
12002) and iBrownl (2004) and we presume th at the crust 
has a completely disordered lattice (lJonesll2004b so that the 
heat transport is contro lled by electron-impurity scattering 
(IItoh&Kohvamalll993h . 

Our neutrino emissivity in the cru st, Ey, in cludes con- 
tributions from bremsstrahlung (Haens el et al. | |1996) and 
Cooper pairing of neutrons in the '5o state dYakovlev et aTl 
[T999>) in the inner crust. Our treatment of neutrino- 
antineutrino bremsstrahlung uses a Coulomb logarithm appro- 
priate for a liquid and does not include band-structure effects 
( iPethick & T horsson 1994; Kaminker et al. 1999), which is 
consistent with our assumption that the lattice is completely 
disordered. In the outer crust we have neglected the plas- 
mon neutrino emissivity ( lltoh et al.lll996h . which is actually 
the dominant mode of neutrino emission there. At the tem- 
peratures of interest, however, the neutrino emission from the 
outer crust is negligible. 

The dominant coolant in the crust is the neutrino emission 
from paired neutrons in the inner crust. For the critical tem- 
perature Tc of the paired neutrons, we choose a Gaussian in 
neutron Fermi momentum kp with a maximum max (^b 7c ) = 
0.8 MeV at kp = 0.8 fm"', and a width 0.28 fm^^ These 
parameters are chosen to reproduc e closely the Tcikp) o f 
lAinsworth et al.l ( ll989l) . As shown bv lCumming et al] (l2006l) . 
the strong temperature sensitivity of the neutrino emissivity 
gy makes the crust temperature rather insensitive to the pre- 
cise depend ence of the critica l temp erature on density. 

Recently, iLeinson & PerezI (12006 ) argued that vertex cor- 
rections in the expansion of the singlet-state interaction would 
reduce Sy by a large fac tor ^ (vp/c)^ <C 1, wh ere vp is the 
neutron Fermi velocity. i Sedrakian et al.l (l2006l) argued that 
although the vertex corrections would change the temperature 
dependence of , the reduction would only be of order unity. 
We use for the rate from lYakovlev et akl (Il999l) . but we will 
include, in §|4] a calculation with this emission suppressed to 
demonstrate how the ignition depth remains sensitive to . 

We couple the steady-state therm al eq uations ([T]i-(l3]) with 
the reaction network calculations (§ 13. lb via an iterative pro- 
cedure. Starting from an initial choice for the temperature 
profile T{r), we integrate the reaction network along the path 
{r[f (r)],p [f (r)]}, where f refers to the Lagrangian time since 
the fluid element was deposited onto the star; in our plane- 
parallel approximation, f(r) = y{r)/m, W\\hy k P / g being the 
column depth. For the convenience of the reader, we show, 
in Figure |2] the column and mass density as a function of 
electron chemical potential pig in the outer crust composed 
of rp-process ashes (§ 13.2.11 ). The curves obey the relation 
jXe °^ (pYe)^^^ °= y'''^, where Ye is the electron abundance. For 
future reference, we also plot in Figure|2]the electron plasma 
frequency {dashed line), which scales with electron chemi- 
cal potential as hcOp^e ~ 78 keV(;Ue/l MeV)^/^. For compari- 
son, the thermal energies are k^T = 43keV(r/5 x 10^ K) < 
hcOp e throughout much of the crust. Our reaction network 



does not account for this large hcOpe in computing (7,n) 
rates. We therefore stop the integration at fie = 20 MeV 
(p > 2 X lO" gcm^^) where the neutron abundance begins 
to rise steeply. The neutrons are not yet in j3 -equilibrium at 
this point, and the density is still less than neutron drip. 
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Fig. 2. — Column density >' = P/g (solid line), mass density p {dotted line), 
and electron plasma frequency fiOp (dashed line) as a function of electron 
chemical potential jl^ in the outer crust where /i^ < 20 MeV. The column y 
is for a neutron star of mass 1.6M., and radius 10.8 km. 

The output of the reaction network gives us the local 
heat Q{r) deposited in the shell {r, r + Ar} per accreted nu- 
cleon. We use Q{r) and the abundance vector Y{r) to solve 
equations ([T])-(l3| and construct a new thermal profile. For 
the inner crust we use the composition and heating from 
iHaensel & Zduiiig (|2003|) ^. We then use the resulting ther- 
mal profile to update the network calculation, and we repeat 
this process until the thermal profile has converged. Typically 
only a few iterations are required, as the heating from electron 
captures is insensitive to the crust temperature. 

3. NUCLEAR REACTIONS IN THE OUTER CRUST 

The evolution of the composition and the associated nu- 
clear energy generation, which produces the flux in the crust 
(eq. 121-1131), is calculated with a nuclear reaction network 
of « 1500 isotopes covering the mass range A = 1-106 and 
extending from proton-rich rp-process ashes to the neutron 
dripline. In this section, we first describe how the electron 
capture and j3 -decay rates are calculated. We next describe 
the heating in the outer crust with two different initial com- 
positions: rp-process and superburst ashes. Building on these 
results, we then develop a method to approximate the heating 
in the outer crust for an arbitrary composition. 

3.1. The Nuclear Model 

We calculate electron capture rates for each time step of the 
integration of the reaction network using a table of electron 
capture transition strengths as a function of excitation energy 
in the daugh ter nucleus and a fast anal ytic phase space ap- 
proximation dBecerril Reyes et al] |2006') that is valid for low 
temperatures. This phase space calculation takes into account 
the temperature and density dependence of the electron cap- 
ture rates and is accurate enough to handle the sharp reaction 

^ The focus of this paper is how the thermal profile, and in particular the 
ignition depth of '-C, changes when the additional nuclear physics is added 
to the outer crust. An exploration of the reactions in the inner crust is left for 
future work. 
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thresholds that make table interpolation impossible. Temper- 
atures and densities are calculated self-consistently from the 
energy generation and composition determined by the reac- 
tion network using an iterative procedure with the thermal 
model discussed in § |2] 

Our nuclear model used to calc ulate the transition s t rength 
functions is described in detail in iMoUer & Randrupl ( Il99(j) 
and here we just give a summary of its properties. We con- 
sider only allowed Gamow-Teller transition strength, which 
is computed in the quasi-particle random-phase approxima- 
tion (QRPA). To obtain the wave functions we solve the 
Schrodinger equation for a deformed folded- Yukawa single- 
particle potential. We then add pairing and Gamow-Teller 
residual interactions treated in the Lipkin-Nogami and ran- 
dom phase approximations. Ground-state deformation param- 
eters a nd masses are obta ined from the finite-range droplet 
model jMoller et al.lll995l) . We calculate neutrino losses for 
each transition, rather than assuming that a fixed fraction of 
{He — Ethi) is lost to neutrinos. As temperatures are low and 
radiative de-excitation timescales are much faster than weak 
interaction timescales, we assume that the parent nuclei are in 
their ground states. 

For the depths considered in this calculation, electron- 
capture-induced neutron emission is not important. We fur- 
ther assume that j3 ^ decay is always completely blocked. This 
is a good assumption in most cases as temperatures are low, 
and because most electron captures either proceed in a double 
step, with the second step proceeding with jXe ^ ^thr.gs-gs. or, 
if they occur at threshold, populate excited states in the final 
nucleus. In both cases, jig is much larger than the decay 
Q-value, which is always equal or less than the ground state to 
ground state electron capture threshold. We will discuss some 
possible exceptions to this in § 13.2.11 

While our main goal is to explore the crust heating within 
the framework of electron capture rates computed from the 
QRPA strength functions, we did adjust the transition ener- 
gies in one case, for which the QRPA predictions are at odds 
with experimental information, and for which this has a strong 
effect on the crustal heating. Our model predicts that the elec- 
tron capture on '""^Rh proceeds to an excited state at 5.3 MeV, 
which woul d result in strong he ating at that depth. Experi- 
mental data jFrevert et al.lll965b shows, however, that while 
^''^Rh is unstable, it also has a small ground-state-to- 
ground-state electron capture branch with a positive Q-value. 
As a result there is no nuclear excitation energy released and 
the heat deposition from this reaction is much smaller than 
predicted from our model. 

In addition to the weak interaction rates, our net- 
work includes neutron capture rates, calculated with 
the statistical Hauser-Feshbach code NON-SMOKER 
(iRauscher & Th ielemann 2000). The corresponding (7,«) 
rates are calculated from detailed balance. For consistency we 
use the NON-SMOKER rates based on masses predicted by 
the finite-range droplet model. For the conditions considered 
here the neutrons are always non-degenerate; denoting the 
degeneracy hy T] — iJ.„{kQT)^^ , we require for degeneracy an 
abundance 

For typical values T = 5 x 10^ K, p = 2 x 10^' gem ^, and 
taking Tj = 0.1 as a fiducial threshold, we have Y„ > 8.4 x 
10^^ as the abundance threshold beyond which neutron de- 
generacy becomes important. At the densities we consider. 



the neutron abundance is well below this threshold. Our the- 
oretical (n, y) and (7,n) rates do not account for the suppres- 
sion of photons below an energy of h(Op_e in a dense plasma. 
Where the neutron separation energies S„ > hcOp e this does 
n ot lead to a large correction in the rate. As we show in 
§ 13.2.11 at jig > 14 MeV some of the nuclei present have suf- 
ficiently low S„ that (7,n) reactions would occur if the pho- 
todisintegration rates were not suppressed. A calculation of 
the (7,n) and (n, 7) rates in this regime is outside the scope of 
this paper; for now, we will indicate their possible impact with 
the recognition that a more careful treatment of these reaction 
rates is warranted. 

The output of the reaction network is the crust composi- 
tion and heat deposition dQ = dEaac + dEg — dEy for each 
timestep dt — {p /m)dr. Here dE = Y.idYi BEi is the nuclear 
energy generation calculated from the abundance change dYt 
and binding energy BEi of each isotope /, dEe — fiedYe is 
the energy change due to changes in electron abundance, and 
dEy is the energy released as neutrinos from electron capture. 
Here all heat depositions are per unit mass. Thermal neu- 
trino losses, i.e., neutrino-pair bremsstrahlung and Cooper- 
pair neutrino emission, are included as cooling terms £v in 
the thermal calculation (cf. § |2]i. We neglect the contribu- 
tion from the change in Coulomb lattice energy dEi^t arising 
from the change in nuclear charge during an electron capture. 
Consider a single ion species (Z, A) with number density riion- 
The lattice energy is mostly due to the Madelung term, E\^( « 
— 0.9(nion/p)^Brr, where F — (Ze)^/ (ak^T) and a is the ion 
sphere radius. The ratio of the lattice energy to electron en- 
ergy is therefore °^ aT?!^ , where a ~ jhc; with the nu- 
merical coefficients inserted this ratio is w 3 x lO^^Z^/^ -C 1 
throughout the outer crust. As a result, the shift in the lo- 
cation of the transitions is small. Moreover, for a two-step 
electron capture, the contribution to dQ comes from the dif- 
ference in li/iiat between the two transitions, i.e., in the rel- 
ative shift of the two transitions. This difference contributes 
w 76 keVu~Hp/10'° gcm-^'fl^A-^l"' < 1 keVu^^ for the 
cases we consider here. 

3.2. Heat Sources in the Outer Crust 

Using the model described above, we computed the heat 
release for different starting compositions. In the following 
discussions, we highlight the dominant transitions in each. 

3.2.1. rp-Process Ashes 

We first calculate the heat sources in the outer crust us- 
ing rp-process ash es from a one-zone X-ray burst model 
dSchatz et al.ll200lT) . The peak of the nuclide distribution is 
around A = 104, and it is useful t o compare our results with 
those of iHaensel & ZduniH (l2003h . who started with a single 
nuclide, '"^Pd. The initial composition is dominated by "^^Cd 
(25% by mass), '"^cd (9.6% by mass), and ''**Ge (7.1% by 
mass). Other mass chains with mass fractions above 2% are 
A ~ 64, 72, 76, 98, 103, and 106. There are also some small 
residual amounts of '^C (0.7% by mass) and '^He (0.5% by 
mass). Figure [3] displays (Z) {lower solid line) and (A) {up- 
per solid line) for this mixture as a function of pie- The rms 
deviations about (Z) and (A) are indicated {shaded bands) 
as well. For comparison, we also plot the Z and A {dashed 
lines) from Haensel & Zdunik (2003). As the fluid element 
chemical potential increases beyond /if = 4 MeV, ^^C fuses 
into heavier elements, with carbon being depleted by 90% at 
pie = 5.8 MeV (p 4.5 x 10^ gcm"^). This leads to the rise 
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in (A) shown in Figure|3] The abundances of '^C is so small, 
however, that they contribute a negligible amount of heating 
at these shallow depths. 
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Fig. 3. — Composition ((Z) and (A); solid lines) for the compre.ssion of 
the rp-process ashes from X-ray bursts. We show the rms deviation about (Z) 
and (A) {shaded bands) for this distribution of nuclides. For comparison we 
also show the Z and A for the composition used bv iHaensel & Zdunik 1 ,200 3„ 
dashed lines). 

We Start our calculation at = 0.38 MeV (p = 6.2 x 
10^ gcm^^). There is an initial spike in heating as the fluid 
element comes into )3 -equilibrium. This heating is an artifact 
of our calculation, but, being located at the outer boundary 
where T is fixed, does not affect the resulting thermal pro- 
file. As the fluid element is moved to greater depth, the rising 
jj.e induces electron captures that reduce (Z) while preserv- 
ing (A) (Fig. IS. For lie > 14 MeV (p > 6 x lO'" gcm^^), 
(A) is no longer constant, as (7,n) and («,7) reactions al- 
ter the relative abundance of different mass chains. Figure]?] 
shows the integrated deposited energy {solid line) as a func- 
tion of /ig. In this plot we set the integrated heat release 
to zero at /ig = 1 .0 MeV so that we can compare our re- 
sults to those of iHaensel & ZdunikI (I2003L dashed line). In 
the range 1 MeV < < 14 MeV, electron captures deposit 
66.3 keVu^', which considerably exceeds the previous esti- 
mate of 16.9 keVu"^ (IHaensel & Zdunikl2003h . In fliis range 
of /ig, neutrinos carry away 37.2 keVu^^ 
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Fig. 4. — Integrated heat release in the crust from the compression of the rp- 
process ashes from X-ray bursts (solid line). For comparison, we also show 
{das hed line) the integrated heat released in the model of Haensel & Zduni^ 
120031) . as well as the results from our approximate model {dotted lines, see 

§113). 



This much larger heat deposition is due to our use of re- 
alistic electron capture rates, which often lead to the popula- 
tion of excited states in the daughter nucleus rather than the 
daughter ground state as assumed previously. Electron cap- 
ture thresholds are therefore increased by these excitation en- 
ergies. More importantly, however, the radiative de-excitation 
of the excited state deposits the excitation energy as heat, 
rather than having it carried off by neutrinos. Note that the 
sum of the heat deposited in the crust and the neutrino loss is 
83.3 keVu^', which is only slightl y larger than the compara- 
ble total heat for the calculation of Haensel & Zdunik ( 2003i) , 
where a single mass A = 106 was used. 

Therefore, electron captures on even-even, odd-even, or 
even-odd nuclei, which tend to occur where the /Xg = isthr, can 
now deposit considerable amounts of heat. These reactions 
were not considered as heat sources before, as it was assumed 
that they proceeded from ground state to ground state. In ad- 
dition, for the electron captures on the odd-odd nuclei — the 
main heat sources in previous work — the fraction of energy 
lost to neutrinos is considerably reduced. 

As an example, a prominent heat source that can be iden- 
tified in Figure]!] is the two step electron capture transition 
of ^""^Ru into "^'*Mo at He ^ 5.2 MeV. The electron capture 
on '"^Ru has a threshold for the transition into the daugh- 
ter ground state of ^thr.gs-gs — 5.09 MeV and is predicted to 
proceed to a low-lying state in ^"'^Tc at E^^c = 0.15 MeV. 
Therefore, this transition occurs at /ig « -Ethr.gs-gs + -Eexc — 
5.24 MeV. The subsequent electron capture on "''*Tc has 
a lower ground state to ground state threshold /sthr.gs-gs = 
1.96 MeV and is dominated by a transition to an excitation 
energy of Eexc — 2.80 MeV in "''*Mo. For the latter transi- 
tion, about 75% of the transition energy /if — -Ethr.gs-gs ^ £^exc = 
0.48 MeV is lost to neutrinos, but the entire excitation energy 
Esxc for both steps is deposited as heat, resulting in a total heat 
deposition of Q = 3.07 MeV per transition. The heat gener- 
ated per accreted nucleon is then QY104 = 8.6 keVu^', with 
F104 — 2.8 X 10^^ being the abundance in the A = 104 chain. 
If the electron capture on '""^Tc proceeded through the ground 
state of '"^Mo, as assumed in previous works, then about 
75% of the entire transition energy /ig — -Ethr.gs-gs — 3.28 MeV 
would be lost to neutrinos and the heat deposition would be 
reduced by a factor of 3.6 to about 0.82 MeV per transition. 
Another prominent heat source in the region jie ^ 14 MeV is 
the two-step electron capture on ^^Fe at /ig « 11 .8 MeV (with 
heat deposition QFgg = 5 keVu^'). 

Typically, transitions in even-A chains, such as the electron 
capture on '""^Ru discussed above, occur as double steps due 
to the odd-even staggering of the electron capture thresholds, 
although exceptions can occur if the transitions populate high 
lying excited states. In contrast, thresholds in odd-A chains 
tend to increase steadily with electron captures proceeding in 
single steps. However, there are also some important two-step 
electron captures in odd mass chains, most importantly the 
transition from "'^Rh into '"^Tc at = 4.4 MeV. While the 
electron capture on "'^Rh has a ground-state-to-ground-state 
threshold of 1 .6 MeV, the transition is predicted to mainly go 
to a 2.9 MeV excited state in '"^Ru. This effectively raises 
the threshold to « 4.5 MeV. Therefore, even though the elec- 
tron capture on ^"^Rh occurs at an jUg near this "effective" 
threshold, the full 2.9 MeV excitation energy is deposited into 
the crust. Following the electron capture onto "'^Rh, jXg is 
then sufficient to overcome the 3.6 MeV threshold (which in- 
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eludes a final state excitation energy of 0.15 MeV) for elec- 
tron capture '"^Ru — + "'^Tc. In total this sequence deposits 
3.8 keVu-' into the crust (Ym = 1.3 x 10"^). 

The mass chains A = 72, 76, 98, and 103 also contribute 
2-3 keVu^^ each in the region /ig < 14 MeV. In total, there 
are about 15 transitions that generate 2 keVu^' or more and 
that are responsible for about 66% of the total deposited heat. 
The remainder of the heating is produced by a larger number 
of less energetic electron captures. At jie }^ 14 MeV, we be- 
gin to see captures of neutrons liberated by (7,«) reactions on 
nuclei with low S„, such as "^^Kr {S„ — 1.9 MeV). As ex- 
plained in § 13.11 the (Y-.n) rates are strongly suppressed where 
S„ < T2(0p,e by the absence of photons with © < (Op,e and so the 
reactions we see are probably an artifact of our network. It is 
interesting to note, however, that in a crust composed of a mix 
of a number of different nuclei with different mass numbers, 
the neutrons released by nuclei with low S„ tend to be recap- 
tured quickly by nuclei with larger neutron capture 2- values. 
It is only at a larger depth, where most A-chains have reached 
low S„, that an appreciable neutron abundance would begin to 
appear. 

Figure |5] (solid lines) demonstrates the effect on the ther- 
mal structure of the crust by the heating from all of the reac- 
tions described above. Plotted are the temperature, in units of 
10^ K (top), and the luminosity L (bottom), which is multi- 
plied by niu/M so that it is in units of MeV. For comparison, 
we also show (dotted lines) the th ermal profiles obtained us - 
ing the composition and heating of Haensel & Zdunikl(l2003h . 
The crust temperature is clearly elevated in the outer crust. 

0.55 F 1 ' 1 ' 1 ' 1 ' 




P/g (g cm") 

Fig. 5. — Temperature (top) and luminosity (bottom) in tlie neutron stai' 
crust, as functions of P/g, for a crust composed of heavy rp-process nu- 
clei (solid lines). For comparison we also show the temperature and lumi- 
nosity (dotted lines) obtained using the crust composition and heating from 
[Haensel & Zdunik (2003). We multiply the luminosity by Htu/M so that the 
plotted quantity is in energy per accreted nucleon. 



An effect that we have not included in our calculations ex- 
plicitly and that could in principle lead to additional cooling 
is an Urea process that might occur for single-step, ground- 
state-to-ground-state electron captures, when the subsequent 
electron capture is blocked until significantly higher density. 
If the temperature is sufficiently high for such an electron 
capture to occur pre-threshold, the reverse )3 ^ decay is not 



blocked. Until /ig has risen to a sufficiently large value the 
electron capture occurs in tandem with a ground-state-to- 
ground-state decay in the reverse direction and an Urea 
process will result. While such transitions are rare (most elec- 
tron captures proceed as double steps, or to excited states) 
there are a few cases where they do occur, particularly in odd- 
mass chains far from neutron sub-shell closures. 

An example is '"^Y, which has a ground-state threshold 
Ethr.gs-gs = 10.83 MeV, while the electron capture daughter 
'"■'Sr has ^thr.gs-gs = 13.39 MeV. For non-zero temperature 
an Urea process can occur for a small range of fie where both 
^"^Y electron capture and "'^Sr decay are fast compared 
to the accretion timescale. Using the electron capture and 
decay rates from our QRPA model we estimate that only a 
few Urea cycles can occur for the peak temperature and ac- 
cretion rate considered here. With each cycle releasing about 
0. 135 MeV per nucleus, we find for this particular case an ad- 
ditional neutrino energy loss of 0.3 keVu^' (using the abun- 
dance in the mass 103 chain F103 ~ 7 x lO^'*). Even if such 
processes occur for about a handful of transitions, the total 
neutrino energy loss from this effect will be negligible. 

The rate of energy loss from neutrinos through this mecha- 
nism is obviously very sensitive to temperature; taking a log- 
arithmic derivative of the rate for the pre-threshold capture 
suggests that the neutrino loss rate is 0= T^. For the specific 
case considered here, the rate is so small that the contribution 
from this process is negligible for T < 8 x 10^ K. It would be 
interesting to explore this process in more detail with accurate 
electron capture and j3 ^ rates for a range of temperatures and 
accretion rates. 

3.2.2. Superburst ashes 

The previous calculation described the evolution of '^C- 
poor ashes produced during an rp-process in an X-ray burst. 
For accreting neutron stars that exhibit sup erbursts, which 
typica lly ignite at a density < 10^ gcm^^ dCumming et al.l 
I2OO6I) . the burning during the superburst sets the initial crust 
composition rather than the X-ray burst. The composition 
in the region between the H/He burning shell and the car- 
bon burning shell is also not well-described by X-ray burst 
ashes, as they do not contain sufficient amounts of '^C to 
ignite superbursts. A better estimate for the initial compo- 
sition in superburst systems might therefore be the ashes of 
stable hydrogen and helium burning, which do contain suffi- 
cient am ounts of carbon at the accretion rates considered here 
dSchatz et al. 200 3b). We therefore perform a second calcu- 
lation, in which we use as an initial composition the ashes of 
steady- state hydrogen and helium burning at the accretion rate 
of O.SMedd dSchatz et alj |2003b). This ash consists mainly of 
52Cr (36% by mass), ^^C (8% by mass), ^^Co+^^Fe (8% by 
mass together), ^**Ni (6% by mass), and ^''Ni (5% by mass). 
We then assume a composition change at iXe = 4.3 MeV 
(p = 1 .7 X 10^ gcm^^) to the superburst ashes of' Schatz et al.l 
(l2003ai) . The superburst ash consists mainly of ''^Ni (35% by 
mass), ^'^Ni (15% by mass), and ^°Fe (14% by mass). 

The electron captures on steady-state rp-process ashes in 
the region 0.45 MeV < jJ-e < 4.3 MeV result in 24 keVu"' 
of heating, while the superburst ashes deposit 32.5 keVu^' in 
the range 4.3 MeV < < 14 MeV. This latter heat deposition 
is slightly greater than half that deposited by the X-ray burst 
ashes in the same range of jig. For jj.^ < 18.5 MeV the energy 
deposition is dominated by just two transitions, ^^Ni ^^Fe 
at He = 9.2 MeV and ^^Fe ^ '^'^Cr at jU^ = 15. 1 MeV. In both 
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these transitions the second electron capture step populates 
relatively high lying excited states at 3.6 and 2.9 MeV respec- 
tively, hence the relatively large heat release. The fact that 
fewer transitions contribute to the heating is not unexpected 
as the composition is more homogeneous than in the case of 
the rp-process ashes. At /i^. > 18.5 MeV there is strong heat- 
ing (34 ke V u^ ^ ) from both neutron-induced reactions and the 
capture ^^Cr ^ ^'*Ti. 

3.2.3. Dependence on composition 

The different heating obtained for rp-process and superburst 
ashes illustrates clearly that the composition matters, when 
one accounts for electron captures into excited states. The 
excitation energies populated by allowed Gamow-Teller tran- 
sitions depend strongly on nuclear structure, in particular on 
shell and sub-shell closures. Figure |6] shows for each par- 
ent nucleus the main excitation energy of the electron cap- 
ture daughter nucleus. Clearly there are huge differences with 
particularly high excitation energies occurring near the shell 
closures. If the initial composition has a large fraction of its 
nuclei that will go through a region of high excitation en- 
ergy, then particularly large heat deposition will occur. Away 
from the shell closures, nuclear deformation tends to frag- 
ment Gamow-Teller strength and suitable transitions are more 
likely to occur at lower excitation energies. 

Although the (7,n) reactions are suppressed where S,, < 
hcOp e, we note that any neutron-induced reaction sequences 
and their energy release will also depend strongly on composi- 
tion. For example, in a mix of nuclei, (7, n) reactions will tend 
to occur at shallower depths in A -chains that reach low neu- 
tron separation energies first (usually only one nucleus exists 
in a mass chain at a particular depth and the heavier, odd-mass 
chains lose neutrons first). In addition to the neutron donors, if 
suitable neutron-accepting nuclei are available in other mass 
chains at the time the neutrons are released, then a net release 
of heat can occur as the composition becomes more neutron 
bound. Depending on the composition, there will be a depth 
where neutron emission would occur without suitable neu- 
tron capture targets available. From there on an appreciable 
neutron abundance would build up. Even if (7,«) reactions 
are completely suppressed such effects might occur once the 
first A-chains have reached neutron drip deeper in the crust. 
In addition, these neutron capture and (7,«) reactions will 
move material to neighboring A-chains with different electron 
capture thresholds and j3 -decay 2-values and thereby induce 
additional weak reactions that otherwise would not have oc- 
curred at that depth. 

3.3. An Approximate Model 

To explore the composition dependence further, and to dis- 
entangle the impact of the various differences to previous 
models more clearly, we present a simplified model of the 
heating in the neutron star crust. In the s pirit of some previous 
studies dHaensel & Zduniklll990ll2003h we limit the model to 
electron captures only (no neutron reactions), we assume that 
electron captures proceed immediately when j^g reaches or ex- 
ceeds the threshold, and we assume that 75% of the electron 
capture transition energy is released in form of neutrinos. In 
contrast to previous work, we do take into account that elec- 
tron captures populate excited states in the daughter nuclei. 
We make the simplifying assumption that only the lowest ly- 
ing daughter state with significant strength is populated. 

To test this model we calculate the heat release as a function 
of depth for our initial rp-process ashes and compare the result 



(dotted line, Figure |4|l with the full network calculations for 
the rp-process burst ashes. Up to jle ~ 14 MeV there is mostly 
good agreement. The slight offset at low is because the tab- 
ulation of the approximate model does not include captures at 
He < 2.0 MeV. The differences at higher jie are due to heat- 
ing from neutron-induced reactions in the reaction network. 
For captures onto nuclei for which jUg ^ fithr our approximate 
model also underestimates the fraction of heat carried off by 
the neutrino. Nevertheless, this approximate model provides 
a good estimate of the crustal heating. 

To explore the composition dependence further. Figure [T] 
maps out the heating, per nucleon, deposited in the crust as a 
function of mass number A at each depth jJ-e, if the crust were 
composed of a single species with that mass number For each 
A, we start with an initial isotope (Z,A) that is stable at the 
low-density boundary. We then use our approximate model 
to trace the heat deposition from successive electron captures 
with increasing jXe- The area of each dot indicates the heat 
deposited in the crust for each transition, and its vertical loca- 
tion indicates the jig at the transition. The largest heat deposit 
is for A = 38 at He = 17.3 MeV with Q = 0.28MeVu-i. For 
this calculation we only include electron captures in the re- 
gion He < 20 MeV. 

In Figure |7] we also show the total deposited heat as a func- 
tion of the mass number A for each mass chain [top panel, 
black histogram). The greatest crustal heating would occur 
for a composition dominated by A = 38 or A = 42 with heat 
deposition reaching 0.45 MeVu^'. Heavy ashes tend to re- 
lease less energy per nucleon as excitation energies are simi- 
lar but the nuclei include more nucleons. Nevertheless, shell 
effects lead to large heating for A = 68-91 with particularly 
high values for A — 70, 84, 86, and 88. In systems without su- 
perbursts, rp-process ashes rich in such nuclei would lead to 
a hotter crust. Superburst ashes in the mass range A = 60-66 
and the rp-process ashes dominated by A = 104 are less favor- 
able for heating the crust. 

For comparison, we also show the heat deposition calcu- 
lated with the same method but assuming ground-state tran- 
sitions only {grey histogram). Taking into account realistic 
transition energies for electron captures leads to a dramatic 
increase in heat deposition that is strongly composition de- 
pendent. 

4. SUPERBURST IGNITION 

As shown in Figure |4] the heat deposited in the outer crust 
is roughly 4 times greater than that computed in an approach 
that neglects captures into excited states. This greater heat 
release in the outer crust raises the temperature (Fig. |5]l in 
the region where '^C would unstably ignite, and hence one 
signature of this enhanced heating is a reduction in the igni- 
tion column depth jign for superbursts. We calculate the col- 
umn depth jign for unstable carbon ignition using the compo- 
sition from steady-state rp-process burning and the superburst 
ashes as described in § 13.2.21 These ashes do contain ^^C at 
mass fraction X('^C) = 0.08. After solving for the thermal 
profile (Fig. [8]l, as described in § |2] we locate the ignition 
column Vi^n usi ng the o n e-zon e stability analysis described 
in Gumming & B ildsten| (l200Jj). Our reaction rate is taken 
from Gaughlan & Fowler ( 198a), with screening according to 
the prescription of iOgata et all (1199 3^. These results over- 
estimate the enhance ment of the reaction rate (Ogata 19971 
lYakovlevet all 120061) . For the case of '^c + fusion and 
the densities of interest, we find that our ignition depths in- 
crease by < 20% when using the screening term indicated by 
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Fig. 6. — For each parent nucleus we show the excitation energy of the dominant state in the daughter nucleus populated by electron capture. Electron capture 
sequences passing through regions with high excitation energies in the daughter will result in increased heat deposition. 



20 24 28 32 36 40 44 48 52 56 60 64 68 72 76 80 84 



92 96 100 104 




20 24 28 32 36 40 44 48 52 56 60 64 68 72 76 80 84 

A 



92 96 100 104 



Fig. 7. — Electron captures in the outer crust. The mass number A is along the horizontal axis; accretion pushes the fluid to greater /i^, which is upward on the 
bottom panel. Each dot on the bottom plot indicates the heat deposited into the crust by a given mass chain at that ^ig, with the area of each dot being proportional 
to the amount of heat deposited. The top panel shows the total heat deposited in the crust (black histogram) for that mass chain, as described in[33] In particular, 
the heat released by (n, 7) reactions is neglected. For comparison, the heat that would be deposited if all transitions were ground-state-to-ground-state is also 
shown (top, grey histogram). 

in the outer crust. The cases shown have the same boundary 
conditions as described in § |2] As Figure |7] shows there is 
a wide variation among the different A-chains in the amount 
of heat deposited and the deposition depth, which leads to a 
sizable difference in temperature at column y ~ 10^^ gcm^^. 

To follow how this variation in temperature affects the ig- 
nition depth, we then computed where a mixture of the the 
crust material and '^C, with XQ^C) — 0.2, would unstably 
ignite. The ignition points are shown in Figure |9](/7/flc^ dots), 
and we plot yjgn as a function of the mass chain A in Fig- 
ure [TO] (^j/acfe histogram). For comparison, we also show the 
ignition depth computed when captures into excited states are 
neglected {grey histogram). There is a general trend that the 
ignition depth decreases with increasing A, even for the case 
with no excitation energy. This is a consequence of two fac- 
tors. First, the electron density at a given column (pressure) 
increases slightly with (Z) to offset the negative Coulomb 
contribution to the pressure. Second, the timescale for a tem- 



lYakovlev et all (l2006l) . 

These ignition points are marked (squares) on Figure[8] The 
ignition column for the four models, from shallowest to deep- 
est, are jign/llO'^ gcm^^) ^ 3 7^ 4 3^ 4 9^ ^j^^j 5 3 ^^it cor- 
responding recurrence times are free = Jign/'w = 5.6, 7.2, 7.4, 
and 9.5 yr. The change in outer crust composition from a 
single mass A —56 with only captures into ground state to a 
realistic mixture of rp-process and superburst ashes shortens 
the recurrence time by about 20%. The increase in temper- 
ature at y < 10^^ gcm^^ is much less dramatic between the 
two crust models. This is due in part because both the steady- 
state rp-process and the superburst produce mass chains that 
are devoid of captures at low pie (cf. Fig^Q. 

We use the approximate model (§ I3.3l l to investigate how 
the ignition depth depends on the heating in, and hence com- 
position of, the outer crust. For each mass chain, we com- 
puted a thermal profile for a crust composed of that nuclide. 
As shown in Figure |9] there is a wide variety in temperatures 
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Fig. 8. — Thermal profiles for crusts with an outer layer of steady state 
rp-process ashes, an outer crust composed of superburst ashes, with heating 
compute d either b y our network calculations (solid lines) or according to the 
model of Haensel & Zdunik 1 1990, dotted lines), and an inner crust computed 
from the model of Haensel & Zdunik 1 1990). The upper set of curves have the 
Cooper-pair neutrino emissivity suppressed. The column where '^C unstably 
ignites is indicated for each solution (squares). 
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Fig. 10. — Carbon ignition depth yig^ for crusts composed of different mass 
chains A. The thermal profile is calculated using our approximate model. The 
difference in heating between the cases for which captures into excited states 
are included (Mack histogram) and neglected (grey histogram) are clearly 
evident. For co mparison, we also show the fitted ignition depths according to 
ICumming et al. 12006, Table 1) (thin dotted lines). 
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Fig. 9. — Temperature in the neutron star crust for mass chains with A = 
21-105. We show two models, one with captures into excited states included 
(dark grey lines) and one with only ground-state captures (light grey lines). 
In each case we set M = 3.0 x lO'^ gs^' and impose the same boundary 
conditions. For comparison, we also show for each model where '~C would 
unstably ignite (dots). 

perature perturbation at a column of 10 gcm^^ to diffuse 
away increases by a factor of « 3 as A increases from 21 to 
105. Both of these effects act to reduce the column needed to 
unstably ignite a ^^C-enriched plasma. This variation with A 
is probably exaggerated somewhat, because we did not self- 
consistently solve the thermal profile with the admixture of 
'^C aty < jign, which would reduce the disparity in the mean 
value of Z and A among mass chains. Nevertheless, this exer- 
cise demonstrates how the superburst ignition depth depends 
on the makeup of the ashes of H and He burning, and the im- 
portance of correctly treating the underlying nuclear physics 
of the crust. The range of ignition columns is more con- 
sistent with the fitted rec urrence times (Fig. [TOl thin dotted 
lines) of G umming et alJ {2006. Table 1) when captures into 
excited states are included, but is still not shallow enough to 
accommodate all of the inferred ignition depths. To reduce 
the ignition column further would require a reduction in the 
co re neutrino em issivity below modified Urea, as suggested 
bv lCumming et al. (2006) . 
If the Cooper-pairing neutrino emissivity were suppressed 



in the inner crust, then the ignition depth is decreased further 
(Fig. fTTT i. As before, the inclusion of electron captures into 
excited states reduces the discrepancy between the calculated 
jign and those fitted from observations. 
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Fig. 1 1 . — Same as Fig. 1101 but for the Cooper-pairing neutrino emissivity 
suppressed. 



5. SUMMARY AND CONCLUSIONS 

We show that allowing for electron captures into excited 
states drastically increases the heat deposited into the outer 
crust of accreting neutron stars. The fraction of the reaction 
energy carried away by neutrinos is sharply curtailed. We 
also find that the amount of heating in the crust is very sen- 
sitive to the initial composition forged during the burning of 
light elements in the neutron star envelope. This sensitivity 
is a consequence of the pronounced shell and sub-shell struc- 
ture of neutron-rich nuclei, which can lead to drastic changes 
in shape, single-particle level structure, and electron capture 
strength distributions, even between nuclei of similar Z, A. 

We have explored how this increased deposition affects 
the ignition depth, and hence the recurrence time, of super- 
bursts. For a realistic composition — steady-state rp-process 
ashes overlying superburst ashes — the inclusion of the excited 
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states decreases the ignition depth by about 20%. For the most 
optimistic case we consider, in which the neutrino emissivity 
from Cooper-pairing is suppressed, the ignition depth for ac- 
cretion at 0.3 Eddington is jign = 3.7 x lO'^ gcm^ , closer to 
the value inferred from observations (jign = 5 x 10^ '-2.7 x 
lO'^ gcm^^; ICumming et alj l2006 l). although a significant 
discrepancy remains unless the core neutrino emissivity is re- 
duced below the modified Urea rate. 

To explore the variation of superburst ignition depths with 
composition, we have constructed an approximate model of 
the heating in the crust. Using this model and including 
the neutrino emissivity from paired neutrons, we find that 
at 0.3 Eddington the superburst ignition depth can be as low 
as 1.3 X 10 gcm^^ if the composition is dominated by nu- 
clei with A = 40-44 or 82-86. The corresponding recurrence 
time for ignition at this depth, with the assumed accretion rate 
and neutron star parameters, is 2 yr. This is roughly consis- 
tent with the observations, without invoking any additional 
physics. We are not aware, however, of a mechanism that 
would, for example, avoid the production of nuclei around 
A « 60 dominating the superburst ashes. 

Our results also indicate that crustal heating is very sen- 
sitive to the underlying nuclear physics. Reliable estimates 
of the structure of very neutron-rich nuclei up to A = 106 are 
needed to constrain the lowest lying electron capture strength, 
which determines the heating from the deexcitation of excited 



daughter states. Nuclear masses are needed to determine re- 
action thresholds and neutron separation energies. Finally, 
masses, jS^decay rates, and reaction rates for the proton- 
and helium-induced reactions on very neutron deficient nu- 
clei during hydrogen and helium burning in the ap- and rp- 
processes are needed to accurately determine the initial com- 
position. More work is also needed to extend our calculations 
into the inner crust to determine the heating there. 
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